Geometry and Regularity of Moving Punctures 



OO 

o 
o 

> ■ 

o^ : 

ON . 

o ■ 
\o ■ 

o ■ 
\o . 
o ■ 

"o '. 
cr. 

u 



Mark Hannam/ Sascha Husa,^ Denis PoUney,'^ Bernd Briigmann,^ and Niall O Murchadha'^ 

^ Theoretical Physics Institute, University of Jena, 07743 Jena, Germany 
Max-Planck-Institut fur Gravitationsphysik, Albert- Einstein- Institut, Am Muhlenberg 1, 14476 Golm, Germany 

'^Physics Department, University College Cork, Ireland 
(Dated: June 16, 2006; July 25, 2007) 

Significant advances in numerical simulations of black-hole binaries have recently been achieved 
using the puncture method. We examine how and why this method works by evolving a single 
black hole. The coordinate singularity and hence the geometry at the puncture are found to change 
during evolution, from representing an asymptotically flat end to being a cylinder. We construct 
an analytic solution for the stationary state of a black hole in spherical symmetry that matches the 
numerical result and demonstrates that the evolution is not dominated by artefacts at the puncture 
but indeed finds the analytical result. 
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Recent breakthroughs in numerical relativity have 
made it possible to simulate the evolution of black-hole 
binaries through several orbits, inspiral, merger and ring- 
down [H, [3, H, 01 ■ It is now possible, after over forty years 
of research, to study black-hole mergers in full general 
relativity, and to calculate the resulting gravitational- 
radiation waveforms. These advances represent a major 
step forward for all of black-hole and gravitational-wave 
astronomy and astrophysics. 

Although these methods have met with spectacular 
success, it is not yet clear how and why they work. In this 
letter we focus on the most popular method, called "mov- 
ing punctures" [1, 0] , and use simple geometrical consid- 
erations to address this question. In the moving-puncture 
method the black holes are conveniently described in the 
initial data in coordinates that do not reach the black 
holes' physical singularities; as the coordinates approach 
each singularity they instead follow a wormhole through 
to another copy of the asymptotically flat exterior space- 
time. These extra copies are compactified so that their 
infinities are represented by single points on a numerical 
grid, which are called "punctures" [1, Q. In this con- 
struction, the wormhole topology is captured by a single 
function, a conformal factor ip, which diverges at each 
puncture. Early "fixed puncture" evolutions 0, fac- 
tored the singularity into an analytically prescribed con- 
formal factor, and the gauge conditions prohibited the 
punctures from moving across the grid. The "moving 
puncture" approach manages to evolve the full conformal 
factor, and the punctures are allowed to move. This 
seemingly minor modification proved to be the last piece 
of the black-hole binary puzzle, and has made long-term 
stable simulations routine for many research groups. 

However, the dynamical behavior of the punctures in 
this method is entirely unknown. Do they continue to 
represent compactified infinities? Does the evolution 
reach a final, stationary state, or do gauge dynamics per- 
sist? And, crucially, does the method accurately describe 
the spacetime, or does it rely on numerical errors near an 



under-resolved puncture, implying that it may fail when 
probed at higher resolutions or for longer evolutions? Fi- 
nally, why do our coordinate conditions accurately re- 
produce angular velocities as would be measured in an 
asymptotic rest frame (see e.g. Q)? 

In this letter, we address these questions in three 
stages. First, we argue that the evolution of the punc- 
tures can be split approximately into advection across 
the hypersurface plus a gauge that locally produces black 
holes in a stationary state. Second, we assume existence 
of a stationary solution and show that a local expan- 
sion for a Schwarzschild puncture matches the numerical 
data. And third, we explicitly construct stationary slices 
for Schwarzschild and discuss their global geometry, to 
which the numerical evolution asymptotes. 

Most significantly, we find that the puncture changes 
character: the slice no longer approaches the other 
asymptotically flat end, but ends on a cylinder of finite 
areal radius. This suggests an elegant new way to repre- 
sent black holes, with asymptotically cylindrical data. 

Let r be the distance to one of the punctures. Internal 
asymptotically flat ends are characterized by 3-metrics 
of the form gij — ip^gij, where (jij is finite and the con- 
formal factor ■0 diverges as ip ^ 1/r. The method of 2| 
introduces a regular conformal factor x — "0"^) while 3| 
uses (j) = logV"- Ignoring the logr singularity, standard 
finite-differencing is used near the punctures. The evolu- 
tion method is based on the Baurngarte-Shapiro-Shibata- 
Nakamura (BSSN) formulation [l^i llH- For the lapse 
we consider "l-|-log" slicing, {dt — (i^dija = —2aK, and 
for the shift the "Gamma-freezing" condition, = 

Puncture evolutions use the standard 3-1-1 decompo- 
sition, and the equations can be brought into the form 
{dt — Cp)u = F, where u is a state vector and F a source 
term independent of the shift. In this context we do not 
have to consider the BSSN variable F* since it is derived 
from cjij. Of special importance is the motion of the 
puncture itself, which is marked by a logr pole in or 
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equivalently a zero in x 1 which are advected by the 
shift Uke all other variables. As noted in 0, H, d, [l^ . 
since the shift can describe arbitrary coordinate motion, 
any shift vector that does not vanish at the puncture can 
generate motion of the puncture. 

The Gamma-freezing shift manages two feats, stabiliz- 
ing the black holes against slice stretching, and moving 
the punctures. Considering that in actual runs those two 
effects appear to happen independently of each other, we 
propose that the shift can be split as — /3* -I- , 
where (31^^ counters slice-stretching and is respon- 
sible for puncture motion. If we have exact stationarity, 
then {dt — -C^^^,, )u = and + F = 0. Note that in 

the early phase of a binary inspiral simulation the vector 
(^)° should be an approximate helical Killing vector, up 
to a rigid rotation. We can choose either corotation, or 
vanishing rotation at infinity, which would imply that the 
punctures' coordinate speeds equal their physical speeds 
as seen from infinity. 

In general, this shift decomposition will hold only ap- 
proximately, but it should be an excellent approxima- 
tion in the immediate vicinity of the punctures. Further- 
more, if the shift is regular at a puncture with a lead- 
ing constant plus higher order terms, = 6q -I- 0{r), 

vection of the puncture without further approximations, 
i.e., u = u(x -\- hot), since CtoU — h^diU. 

We conclude that, assuming approximate stationarity, 
the question of numerical stability of the moving punc- 
ture process can be approached by splitting it into a 
standard advection problem and a stability analysis of 
a stationary puncture solution. The main open issue is 
whether there exists a regular, stationary solution for a 
single puncture that the method can find. Since the key 
novel features of this stationary solution already occur 
for a non-moving, spherically symmetric puncture (com- 
puted with the moving puncture method), we focus on 
this case here and leave the general case to future work. 

Consider initial data for a single puncture with mass M 
at i = 0, with a = 1, = 0, gij = A.j ^ 0, K ^ 0, 
and 00 = log(l + When inserted into the BSSN and 
gauge equations, the data evolves and develops certain 
powers of r at the puncture [13] . If a regular, stationary 
state is reached, then all variables should possess power 
series expansions at r = that satisfy Cp^i^u -I- F = 0. 
We make the following ansatz for the single, non-moving 
puncture case in spherical symmetry: 
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Here r = {x^ + + z^)^/^ is the coordinate radius of 
quasi-isotropic Cartesian coordinates (x, y, z). Note that 
r is continuous but not differentiable, the radial nor- 
mal vector = x'^/r is discontinuous at r = 0, and 
diO{r°) — 0{r^^). The evolution equations and the con- 
straints result in 8 independent equations for the 10 coef- 
ficients {pi,p2, Aq, Ai, Kq , i^i , ao , ai , 6i , 62 ) , with M and 
7] free parameters. The Gamma-freezing shift condition 
does not give an additional condition for a stationary so- 
lution. A first result is that ai and Kq ^ are 
required for nontrivial stationarity; a more regular solu- 
tion is not consistent with the equations. For consistency 
ao = 0, that is the lapse has collapsed at the puncture. 
Further simple relations are hi = 2Kq and &2 = —SAqQi 
from the evolution equations for the lapse and the metric. 

Therefore, assuming existence of a stationary solution 
with some minimal regularity at the puncture implies 
specific predictions for kinks and discontinuities in our 
variables. We now compare the ansatz to numerical re- 
sults for a Schwarzschild puncture evolved tot = 50M, at 
which time the data appears to have reached a stationary 
state to within a few percent in all variables. The runs 
were performed with the BAM code, using fixed-mesh re- 
finement and fourth-order accurate finite-difference sten- 
cils [1^. For a central resolution of M/128 the run is 
well within the convergent regime. When zooming into 
the data at this resolution. Fig. [1] we do find the kinks 
and discontinuities consistent with the ansatz. Although 
closer inspection shows some numerical artefacts near the 
puncture (in particular in Aij and the derivatives), they 
remain localized and comparatively small, and the nu- 
merical data are approximated well by the ansatz. 

The most intriguing implication of the power series 
expansion is that the singularity in the conformal factor 
changes during evolution. With oq = 0, 

dt^ = /3'9,0 + ^d,(3' + 0(r) ~ birdr^ -f y . (7) 

For the conformal factor of the initial data, we have 
rdr(f>Q — —1, and dt(f>Q — — &i/2. With 61 7^ we must 
conclude that the evolution of the conformal factor will 
stop only when rdr4> — ~l/2- In other words, 

^0 = O(l/r) evolves into ip^Oil/Vr), (8) 

which is built into ([1]) . The puncture still marks a coordi- 
nate singularity, but since the areal radius of the sphere 
r — is Ro = limr_>o ^"0^ = the slice no longer 

reaches the other asymptotically flat end of the Brill- 
Lindquist wormhole, but ends at a finite areal radius, 
Rq ~ 1.3M, with pi determined numerically. Figure [2] 
illustrates the evolution of a Schwarzschild black hole. 
We have plotted Schwarzschild R versus proper distance 
from the horizon. The initial slice connects two asymp- 
totically flat regions, but during a l-l-log/F-driver evo- 
lution the numerical slice loses contact with the second 
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FIG. 1: Numerical evolutions of a spherically symmetric 
puncture result in specific kinks at the puncture (left). The 
figure on the right represents a true numerical experiment. 
We compare two spherical slices by plotting a as a function 
of K. The dots are the data from the 3D numerical evolution, 
the solid line was obtained by independently integrating up 
Eq. (9). At the puncture, a ki 0.0 and K « 0.3. 




flat, i.e., R ~ I, must pass through such a singular point. 
The only way of resolving this difficulty is for the nu- 
merator of (O to simultaneously vanish at the "singu- 
lar" point. We have two options: either (a) R' = 0, or 
(b) 3M/R - 2 + 2R'^ = 0. Let us discuss them sep- 
arately. In the case (a) we have R' — 0, R ~ 2M. 
Therefore the slice passes through the bifurcation sphere. 
This is the standard moment of time symmetry slice 
through the Schwarzschild solution, which obviously sat- 
isfies the stationary 1-l-log equation because /3 = and 
K = 0. In case (b), the case we are really interested 
in, we can solve the pair of simultaneous equations to 
give R' = VlO - 3, i? « 1.54M. This solution corre- 
sponds to two slices in the Schwarzschild solution. These 
are mirror images of each other, one in the upper half 
plane, one in the lower. These slices do not continue into 
the left quadrant of the extended Schwarzschild solution. 
Rather they asymptote to a cylinder in the upper (lower) 
quadrant of fixed radius. This agrees with the numerical 
observation that the singularity in the conformal factor 
changes from 1/r to \l\pr. These three are the only 
asymptotically regular solutions of the stationary 1-1- log 
equation, up to isometrics. It is easy to show that the 
actual slice exponentially approaches a cylinder of radius 
i?o, i.e., i? « i?o + AexpBZ as I — > — oo, where 

B = 2Rq-^ (3Af - 2i?o) (2M - R^Y^ . 



FIG. 2: Schwarzschild coordinate R versus proper distance 
from the (outer) horizon. The left panel shows the slices at 
t = 0, 1,2,3M, and the right panel shows the slice a,t t — 
50M. The final numerical slice terminates at i? « 1.3A'/. The 
vertical line indicates the horizon at R — 2M, and the six 
points represent x/AI — 1/40, 1/20, 1/8, 2, 5, 8 on each slice. 



asymptotically flat end. It eventually reaches a station- 
ary state and terminates at i? « 1.3Af . 

We now show by explicit construction in the proper dis- 
tance gauge (which is regular through the horizon), that 
such stationary slices exist globally. Stationary l-flog 
slices of Schwarzschild are also considered in a differ- 
ent context in [l6| . The proper distance coordinate is 
denoted as I, dif — /', a is the Killing lapse, /3 the 
Killing radial shift component and R the areal radius (the 
Schwarzschild radial coordinate). On any spherical slice 
through the Schwarzschild solution one finds R' — a and 
— pi^ = \ — 2M/R, and the equation for the station- 
ary H-log slices of Schwarzschild, /3a' = 2aK, becomes 
(using, for example, fv^ ) 



R" 



2R' 



3M _ 9 



2R'- 



i?'2 - 2R' 



(9) 



It is clear that the right-hand side of ([9]) is singular when- 
ever 2M/R - 1 + - 2R' = 0. It can also be shown 
that any solution of ^ that is suitably asymptotically 



The value of the lapse at the horizon, which can be used 
as a simple horizon-finding or merger-time criterion, eval- 
uates to a{R = 2M) « 0.376. It is also possible to pro- 
duce an algebraic solution in terms of an implicit equa- 
tion for a and R; determining i?o ~ 1.31241M from 



3sinh-^(3) = log 



128(2- 

^ M' 



10-3. 



For comparison we have also studied harmonic and 
maximal slicing. With harmonic slicing, (which can be 
generalized to the Kerr spacetime 18|, ll9|), one again 
finds a moment of time symmetry slice and two mirror 
copies of a slice that puncture evolutions driven toward 
stationarity should approach. The slices hit the singular- 
ity a,t R = 0, K blows up as i? ^ 0, and at the horizon 
a{R = 2M) = 1/2. Fixed-puncture evolutions have used 
1+log slicing without shift term. Stationary slices must 
then be maximal 



21|, 123, of the "odd lapse" type, 
a2 ^ 1 - 2M _^ C^^ For C = 3V3MV4 such shces will 
again approach a cylinder, of constant Rq ~ 3A//2. 

Discussion. — Moving-puncture evolutions of a 
Schwarzschild black hole approach a stationary slice that 
neither reaches an internal asymptotically flat end nor 
hits the physical singularity, as might be expected for a 
stationary slice with non-negative lapse Rather, the 
slice ends at a throat at finite Schwarzschild radius, but 
infinite proper distance from the apparent horizon. This 
changes the singularity structure of the "puncture" . It 
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is still a puncture in that there is a coordinate singular- 
ity at a single point in the numerical coordinates, but it 
does not correspond to an asymptotically flat end. In the 
course of Schwarzschild evolutions we have found that the 
throat does collapse to the origin. Where one would have 
expected an inner and an outer horizon, we find only one 
zero in the norm of (^)°, corresponding to the outer 
horizon. An under-resolved region does develop in the 
spacetime (it is the region between the throat and the in- 
terior spacelike infinity) , but we are pushed out of causal 
contact with it. The throat itself has receded to infinite 
proper distance from the outer horizon. Matter fields or 
gravitational radiation will be trapped between the inner 
horizon and the throat, because unlike the gauge their 
propagation is limited by the speed of light; this issue is 
left to future work. 

The main result about regularity is contained in the 
specific kinks and discontinuities discovered in both the 
power series and the numerical data, which are implied by 
our stationary solution. The construction of an explicit 
power series solution in r for a single Schwarzschild punc- 
ture can be repeated for a single puncture with Bowen- 
York linear momentum. Preliminary results indicate that 
the shift can acquire the constant leading term respon- 
sible for the motion of the punctures, and that there 
are no new fundamental regularity issues. The power 
series ansatz matches nicely the numerical data for a 
Schwarzschild puncture. The fourth-order schemes (both 
centered and upwind) in use have been seen to be reason- 
ably successful in this context, the lack of differentiabil- 
ity not withstanding [3, Q • Incidentally, the introduction 
of X results in O(r^) terms and not 0(r'*) as expected, 
but this still results in significantly cleaner data than the 
logr method. Perhaps most importantly, our analysis 
suggests a concrete remedy if the remaining numerical 
issues at the puncture create problems in simulations of 
black-hole binaries, namely to resort to finite differencing 
that is expertly adapted to the discontinuities at hand. 

Our results suggest several directions for future re- 
search directly relevant to the black-hole binary problem, 
such as perturbations of the stationary solutions (includ- 
ing constraint violating perturbations to check constraint 
stability of evolution systems), the clarification of numer- 
ical issues at the discontinuities, and the construction of 
initial data adapted to stationarity, e.g. of asymptotically 
cylindrical data. In the Schwarzschild case local proper- 
ties of the stationary solution allow one to directly read 
off spacetime properties from a numerical solution, e.g., 
the puncture value of K determines the mass; an exten- 



sion to two moving, spinning black holes would be very 
valuable in numerical evolutions. 
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